Hybridization wave as the 'Hidden Order' in URu2Si2 
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A phenomenological model for the 'hidden order' transition in the heavy Fermion material 
URu2Si2 is introduced. The 'hidden order' is identified as an incommensurate, momentum-carrying 
hybridization between the light hole band and the heavy electron band. This modulated hybridiza- 
tion appears after a Fano hybridization at higher temperatures takes place. We focus on the hy- 
bridization wave as the order parameter in URu2Si2 and possibly other materials with similar band 
structures. The model is qualitatively consistent with numerous experimental results obtained from 
e.g. neutron scattering and scanning tunneling microscopy. Specifically, we find a gap-like feature 
in the density of states and the appearance of features at an incommensurate vector Q* ~ 0.67r/ao. 
Finally, the model allows us to make various predictions which are amenable to current experiments. 



PACS numbers: 

The heavy Fermion material URu 2 Si2 exhibits an 
enigmatic 'hidden order' (HO) transition at a temper- 
ature T HO = 17.5K Although being more than 
two decades old, the nature of the transition (which is 
a typical mean-field second order transition) and the or- 
der parameter which describe it are still unknown 0, Hj] ■ 
Various theories have been put forward to explain the 
transition, some of them describing it as a localized one 
111 and some describe it as an itinerant phenomenon 
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In the last few years various new experimental tech- 
niques, previously unavailable, were presented, and now 
provide additional valuable hints as to the origin of the 
Hidden Order transition. The key experimental observa- 
tions relevant for our discussion are: 

(i) In inelastic neutron-scattering measurements fl8j 
(as well as other measurements (l9||) enhanced scattering 
is observed at the incommensurate wave vectors Q* ~ 
0.6,1.4 n/ao below the HO transition. A spin-gap-likc 
feature of A ~ hmeV is observed, as extracted from the 
heat capacity measurements. In many heavy-Fermion 
systems the spin-gap is indicative of the accompanying 
charge gap [201 ] . 

(ii) Angle-resolved photoemission spectroscopy mea- 
surements indicate that the heavy /-band is very weakly 
dispersive, crosses the Fermi energy at the HO transition, 
and settles slightly below it plj . 

(iii) A Fano line-shape in the density of states (DOS) 
which starts to form at ~ 120K has been seen in recent 
scanning tunneling microscopy (STM) measurements [22I 
[23( | . This Fano line-shape is fully developed at ~ 20K 
(inset of Fig.^a)). 

(iv) At the HO transition temperature, the bottom 
of the Fano line-shape develops a gap-like feature which 
evolves with temperature (22 . |23| (inset of Fig. |2£b)). 



Both the Fano parameters and the gap structure depend 
on the STM tip position and modulate in space with the 
positions of the different atoms. 

(v) The hole band develops a hybridization feature 
below the HO transition, corresponding to momentum 
Q w 0.3 it Jan, close to where the heavy Fermion band 
crosses it [22j . 

These observations seem to suggest that a descrip- 
tion of the HO in terms of reorganization of the local- 
ized degrees of freedom may be incomplete. Guided by 
the experimental findings, we construct a phenomeno- 
logical theory for the HO transition. A one-dimensional 
model system is taken for simplicity, which contains a 
light hole (d) band and a heavy electron (/-electrons, 
dashed line) band, depicted in the inset of Fig. [TJ We 
also performed the same calculation using bands gener- 
ated by local density approximation (LDA) as input (to 
be shown elsewhere). Both from the model and from 
LDA based bands we take the bands crossing at momen- 
tum points Q pa ±0.3 ( along (1,0), (0,1) in units of 
tt/uq, where ao is the lattice spacing, in which all mo- 
menta are measured hereafter). The holes first hybridize 
with the local part of the /-electrons, giving rise to a 
Fano line-shape in the DOS [HI, [24-27 1 at temperatures 
that are much higher then the HO transition. The band- 
structure gives rise to enhanced hybridization between 
electrons with momentum Q and holes with momentum 
— Q (and vice versa, Fig. [T|) at specific points along the 
Fermi surface ("hot spots"). The resulting electron- hole 
coherence (which may be considered as an indirect exci- 
ton) is the HO order parameter. Since the "hot spots" 
are located at the band crossing, this incommensurate 
scattering generates a gap structure in the DOS with a 
line-shape similar to that observed in STM experiments 

The starting point of the model is a two-band hamil- 
tonian, 
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FIG. 1: Schematic representation of the band structure as- 
sumed in the model and the incommensurate scattering. The 
Heavy electron /-band and the light hole d-band cross at 
Q m ±0.3 (inset). The resulting scattering hybridizes elec- 
trons with momentum Q and holes with momentum — Q (and 
vice versa), generating a hybridization wave which is the HO 
order parameter. The Q* = 2Q ~ 0.6 momentum transfer is 
in accord with Neutron scattering experiments [l8l ]. 



where ct (/^) creates a hole (electron) with momentum 
k. For simplicity (and as it plays no important role in 
the physics we describe) the electronic spin is omitted. 
In real materials the d — / hybridization sets in at much 
higher temperatures then HO transition and hence the 
remnant heavy band and light hole bands we use here 
are already hybridized. We assume that lighter band is 
more d-character and heavier band has more /-character 
near the T point where we concentrate. 

The origin of the Fano lineshape has been a subject 
of recent study 



- 281] . In the absence of cither trans- 
lation invariance breaking or level broadening, a Fano 
lineshape should not appear [29| . A simple mechanism 
for its appearance is based on local scattering between 
the d-electrons and the /-electrons under the STM tip 
(the effect of the tip may not be negligible, [I(J). The 
hamiltonian is then given by tip = Vq c\.f r= o + h.c. , 
where f r= o creates an /-electron on the lattice site r = 0, 
i.e. under the tip. Alternatively, a similar result as below 
can be obtained by assuming full decoherence between 
the /-orbitals which lie on different lattice sites (29j . 

It is now straight forward to write the Dyson's equa- 
tions, with the solutions for the Green's functions for the 
/- and d-electrons / - 
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- -V^ Xo fffc '9k'' ■ -Here g k is the bare d-electron prop- 
agator, and its susceptibility \o = J2k 9k = — To(i — q), 
which defines Tq (proportional to the bare <i-band DOS, 
which we assume constant for simplicity, and set it as 
To = 1 hereafter ) and q (which describes the ratio be- 
tween the real and imaginary part of the bare self-energy, 
and is known as the Fano factor). The level £o is not the 
bare /-level, but rather a renormalized value, which may 
arise due e.g. to hybridization with other bands, inelas- 
tic scattering effects, interaction rcnormalization etc. Eq 
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FIG. 2: (a) Local DOS for a d-band coupled to the local- 
ized /-electrons, exhibiting an asymmetric Fano line-Shape 
(solid line, see text for parameters). The dashed line is the 
local DOS in the presence of the hidden order. These are in 
good agreement with the DOS (proportional to conductance) 
seen in STM experiments, shown in the inset (extracted with 
permission from [12]). (b) Change in DOS due to the HO 
coupling, Sp = p(V) - p(V = 0), for V = 0.2,0.4, 1. The 
main features are the appearance of a gap structure (which 
scales with V, see Fig.[3]), and a peak structure within the gap, 
appearing for large enough V (corresponding to low enough 
temperature). Both these features have been observed exper- 
imentally [23 . [231 . as seen in the inset (extracted with per- 
mission from 1231). 



can be lower than the original /-level [28(, and there- 
fore, a negative value of eo does not imply an assump- 
tion that the /-level lies below the Fermi energy. The 
local DOS is given by p = — k , gk.k' ■ To be spe- 

cific in comparing our results with experiments we fix 
parameters used as Ti = VqTq = 5 (all energy scales 
are taken in meV), the bottom of the /-band Eq = — 5 
(shifting this energy does not alters the line-shape) and 
the Fano factor q = 1.5 (these numerical parameters are 
maintained hereafter unless stated otherwise). This line- 
shape is to be compared with that obtained from STM 
measurements, shown in the inset of Fig. [2] [22|, [23| ■ In our 
treatment, temperature only enters via the temperature- 
dependence of the couplings. The hybridization feature is 
of A ~ 5mcV~ 6QK, and the HO transition is at 17.5K, 
therefore the finite temperature should smear the DOS 
line-shape, but not change its features. 

The next step is to define the effective interaction be- 
tween the /-electrons and the d-holes which gives rise to 
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the HO. The central part of our proposal for the HO is 
to assume that the most important electron correlations 
affect matrix element between d-band holes with momen- 
tum —Q and f -electrons with momentum Q (and vice 
versa), with Q ps 0.3 (the position where the bands cross, 
Fig- [J). This assumption is consistent with the observed 
sharp onset of the hybridization feature at Q* = 2Q seen 
in quasiparticle interference 



221 ] . Thus, we take the in- 



teraction Hamiltonian as 



rlj =Y,Ukk'ci^ Q fi_ Q f k , +Q c k+Q + (Q ^ -Q) , (2) 



where Uk is strongly peaked around k = (the width 
of Uk determines the range of the effective d — f in- 
teractions). To continue analytically, wc set ft = k' 
and Uk ~ U5k,o, resulting in the effective Hamiltonian 
Hi ~ U Q -Qcl Q fl Q f Q CQ + (Q -4 -Q) . The physical 
meaning of these assumptions is that (i) the most impor- 
tant scattering is with momentum transfer O (l7, 18. 31 [. 
and (ii) the electron interactions are long-ranged [32 1. 
Since the HO transition seems like a mean-field transi- 
tion (based on the shape of the heat capacity @ and the 
dependence of the gap on temperature [23j ) , it is natural 
to decouple Hi n a mean-field way. We thus identify the 
HO order parameter V as 



V = V HO = U- Q , Q (cl Q f Q ) 



(3) 



The 'hidden' nature of this excitonic order is evident from 
the observation that the expectation value of the spin 
and of the charge (density waves) would be zero in this 
state hence there is no primary single spin and or charge 
order in this hybridization wave state. We do expect a 
charge density modulation at momentum Q* = 2Q as a 
secondary effect, which can be however very small (33| . 
In the mean-field approximation, V dependens on tem- 
perature as V oc (Tho — T) 1 / 2 . The interaction hamil- 
tonian in the mean-field approximation takes the form 
Uimf = Vd_ Q fQ + h.c. 

We now write the Dyson Equation for the Green's func- 
tion in the presence of Hi. We point that if Hi would op- 
erate on the bare hamiltonian, then no correction (other 
than at the points k = ±Q) would be observed. However, 
the d-band is already dressed by the Fano interaction de- 
scribed above, and thus the correction is to all momenta. 
Since the bands are symmetric in k, we can treat only 
the — Q — > Q (the full calculation which includes also 
Q — > —Q scattering gives similar results). The Dyson 
equations are 



Gk,k' = 9k, k 

Fq,q = 



■' + V 2 9k-QFQ t QG-Q,k 

v 2 Iq,qG-q-q f q,q 



(4) 



which are easily solved by simple algebra. 

The resulting local DOS p — J^k k' Gk,w is plotted 
in FigHKa) (dashed line) for V = 0.8 and e Q = -2.5. 



in Fig. [2jb) we plot the change in DOS due to the HO 
coupling, Sp = p(V) - p(V = 0), for V = 0.2, 0.4, 1. 
As seen, a gap-like feature develops, its width and posi- 
tion depending on the value of V (which in turn depends 
on temperature). Note also the additional peak at the 
bottom of the gap, which was observed experimentally 
[22T I23I ]. For comparison, in the inset of Fig. [2jb) we 
plot the experimental DOS taken from [22| ■ The qualita- 
tive agreement between the data and theory is evident. 
A quantitative fit would require additional parameters to 
account for the complicated background DOS seen in the 
experiments. 

From the DOS one can extract the size of the energy 
gap A as a function of V (for instance, by taking the 
distance between the peak just left of the gap and the 
point where the curvature is maximal just right of the 
gap). The dependence of A on V 2 cx (Tho — T) is plotted 
in Fig. [Hta). We find that A oc (T HO - T) v , with v sa 
0.4 (solid line in Fig. [HJa)), in good agreement with the 
exponent extracted from the experimental data [23[ and 
in agreement with the mean-field treatment assumed. 

A pronounced feature of the DOS is the appearance of 
an additional peak structure within the gap in the DOS 
(Fig. [21(b)). In the experiment 22, 23| this feature only 
appears at T « Tho which we interpret as large enough 
V, as it grows below Tho- In Fig. Elb) we have plot- 
ted the position of the minimal point of the DOS (E m i n , 
empty circles) and the position of the peak within the 
gap (-©peak, full circles) as a function of the HO order pa- 
rameter V (which, in experiment, corresponds to temper- 
ature) . The peak structure only appears for large enough 
V, and scales with temperature with the same exponent 
as the energy gap itself, E pca ^ oc A oc (Tho — T) 0A . 
-E m i n , on the other hand, quickly saturates to become 
temperature-independent. Both these quantities could 
be examined in a detailed STM experiment. 

In the inset of Fig. [3^b) we have plotted the DOS 
change 8p for different values of the Fano factor, q = 
0.5, 1, 1.5, 2 (for V = 0.5). In the STM experiments, the 
value of the Fano factor changes, depending on whether 
the tip is positioned above a Si site or a U site. Presum- 
ably, when the tip is above the Si site it has better cou- 
pling to the d-band, which effectively increases the Fano 
factor. Regardless of the reason for the g-factor modu- 
lation, we predict that the height of the peak structure 
within the gap will also modulate (and will be commen- 
surate with the modulation of q) , but its position will not 
change (dashed line in the inset if Fig. [3](b)). 

The order parameter we discuss here Vho = 
U-q^q (c^qfo^ also carries a nonzero expectation value 
for the dipole moment modulation at Q*. In the coher- 
ent excitonic state |\E r ) = rL( M fc c -fc + |0) there is a 
nonzero expectation value for the dipolar moment at mo- 
mentum Q*: (|-P2q|) = ('3 / |re l2 '2''|iI'), as long as the par- 
ity of d-&nd /-orbitals is different. This modulation of the 
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FIG. 3: (a) The energy gap, as obtained from the DOS 
(Fig. 0) , as a function of the HO order parameter V squared, 
which is proportional to (Tho —T). The solid line is the func- 
tion A oc (Tho — T) 0,4 . Inset: the dependence of Tho on 
the interaction strength U, indicating a possible first order 
transition at U = U c . (b) Position of the minimal point of 
the DOS (-Emm, empty circles) and the position of the peak 
within the gap (i? pGa k, full circles) as a function of the HO or- 



der parameter V' 



(which corresponds to temperature). 

0.4 



While 

i?min saturates at large V, -E pca k ~ V ~ |T — Tho| u ' 4 in the 
same way as the gap A (Fig. |3J. Inset: change in the DOS, 
8p for various values of the Fano factor q (for V = 0.5). 



dipolc moment could be observed with probes that cou- 
ple to electric density modulation, for example an STM 
with a ferroelectric (or hybrid ferroelectric-metallic) tip. 

The formalism described above allows us to calculate 
the critical temperature of the HO transition. Working in 
a Nambu-like space mixing d- and /-electrons and using 
Eq. [3J we identify the self-consistency equation for the 
HO, 



1 



U- 



= T H0 ^2fQ{iu>n,Q)g- 



)(iw„) 



(5) 



where fQ(iui n ) and g-Q(iuj n ) are the Matsubara Green's 
function in the presence of the Fano interaction, evalu- 
ated at momentum Q, and ioj n are Matsubara frequen- 
cies. We have calculated the critical temperature Tho 
for the parameters above as a function of U, shown in 
the inset of Fig. Ela). We find that for U below a certain 
value U c , Eq. [5]docs not have a solution, due to the com- 
pact nature of the interaction in momentum space. This 



feature may explain why the hidden order has not been 
universally observed in other heavy-Fcrmion materials, 
even ones with a similar band structure. It may also 
be relevant to the destruction of the HO and appear- 
ance of magnetic ordering with applied pressure, with 
a 1st order transition between them [34, 35 1 A possible 
scenario for the connection between the HO and mag- 
netic order (clearly avoided in the present study which 
discusses spinless electrons) is the following. As pres- 
sure is applied the interaction strength decreases, un- 
til finally it becomes lower than U c . The magnetic in- 
teraction between the localized /-electrons, previously 
screened by the conduction electron which are coupled 
to the /-electrons in the HO state, becomes significant 
and gives rise to magnetic ordering. This scenario will 
be studied in detail in the future. 
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